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In  recent  years,  great  advances  have  been  made  in  the  area  of 
computerized  structural  analysis  and  optimal  design.  One  class  of 
optimization  techniques  called  mathematical  programming  methods 
have  been  applied  extensively  to  the  design  of  structural  systems. 
In  most  cases  the  resulting  problems  are  nonlinear  in  the  design 
variables  and  thus  nonlinear  programming  techniques  are  required 
for  the  solution.  The  difficulty  with  these  methods  is  that  as 
the  number  of  design  variables  becomes  large  . . on  the  order  of  50 
or  more  ..  problems  with  computer  costs  (or  time)  and  numerical 
accuracy  become  significant,  thus  limiting  their  applicability. 

The  goal  of  the  work  presented  here  was  to  develop  an  efficient 
optimization  technique  for  the  design  of  large  structures. 

While  nonlinear  programming  methods  are  impractical  for  large  prob- 
lems, linear  techniques,  if  applicable,  can  still  be  used.  This, 
then,  was  the  general  approach  taken  in  the  current  research. 
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2. 


The  use  of  linear  programming  methods  in  the  design  of  trusses 
has  been  investigated  by  a number  of  researchers.  Dorn,  Gomory,  and 
Greenberg  [5]  have  shown  that  the  optimum  design  of  a truss,  sub- 
jected to  stress  limitations,  can  be  formulated  as  a linear  program- 
ming problem.  The  authors  considered  the  geometry  of  the  structure 
to  be  fixed  and  took  the  member  areas  as  design  variables.  All 
trusses  obtained  from  their  approach  were  determinate,  since  only 
conditions  of  equilibrium  were  imposed.  Farshi  and  Schmit  [6] 
developed  a technique  for  the  optimization  of  trusses  based  on  the 
force  method  of  analysis.  Initially  only  stress  constraints  and 
the  conditions  of  equilibrium  were  considered,  and  then  in  an  itera- 
tive manner  the  conditions  of  compatibility  were  added  until  a 
suitable  design  was  obtained.  In  studies  by  Reinschmidt, 
Norabhoompipat  and  Russel  [12,13],  an  approach  similar  to  that  of 
[6]  was  used.  One  conclusion  from  their  work  was  that  the  approx- 
imate solution  found  by  satisfying  only  the  equilibrium  conditions 
was  close  to  the  optimum  design  obtained  by  enforcing  the  condi- 
tions of  compatibility  as  well.  In  an  investigation  by  Schmit  and 
Farshi  [16]  an  efficient  technique  for  the  optimization  of  trusses 
subjected  to  both  stress  and  displacement  constraints  was  developed. 
In  order  to  obtain  the  optimum  design  a sequence  of  linear  program- 
ming problems  were  solved.  The  optimization  algorithm  was  based 
on  the  method  of  inscribed  hyper spheres  [1] . High  efficiency 
was  achieved  by  using  several  approximation  concepts  such  as 
design  variable  linking,  temporary  deletion  of  noncritical  con- 
straints, and  Taylor  series  expansions  for  response  variables 
in  terms  of  the  design  variables.  The  authors  concluded 
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from  this  study  that  the  use  of  approximation  concepts  along  with  j 

i 

■* ' linear  programming  techniques  could  lead  to  a new  generation  of 

efficient  and  practical  optimization  procedures. 

• • 

The  design  of  frames  has  also  been  formulated  such  that 
^ linear  programming  techniques  can  be  applied.  Moses  [11]  performed 

successive  Taylor  series  expansions  to  obtain  linear  design  con- 
straints. An  iterative  approach,  based  on  the  simplex  method,  was 
used  in  solving  for  the  final  design.  In  the  work  of  Rubinstein 
and  Karagozian  [14],  plastic  hinge  theory  was  used  to  generate  the 
equations  of  constraint  on  the  moment  carrying  capacity  of  the 
members.  The  design  criteria  was  safety  against  collapse  and  a 
limitation  on  lateral  deflections  in  the  elastic  range.  Romstad 
and  Wang  [15]  used  another  approach  to  formulate  the  optimization 
of  frames  as  a linear  programming  problem.  Linear  behavior  con- 
straints on  stresses  and  displacements  were  developed  by  initiating 
local  changes  in  the  design  parameters  and  then  determining  the 
resulting  stress  and  displacement  redistributions.  A design  cycle 
in  this  formulation  involved  analyzing  the  structure  using  assumed 
design  variables  via  the  displacement  method  of  structural  analysis. 

These  results  were  then  used  to  generate  the  linear  behavior  con- 
^ straints  which  comprised  the  linear  programming  problem.  The  solu- 

i . tion  of  this  linear  programming  problem  resulted  in  a lower  weight 

I design  with  modified  design  variables,  which  completed  the  cycle. 

In  the  work  presented  here,  the  optimum  frame  design  problem  is 
formulated  in  a new  way  such  that  linear  programming  methods  can 
more  efficiently  be  applied.  Stress  constraints,  as  well  as  conditions 
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of  equilibrivun  and  compatibility  are  considered.  The  approach 
does  not  require  a separate  time  consuming  analysis  phase  to  deter- 
mine member  stresses  or  forces  (as  does  the  Romstad's  and  Wang's  | 

i 

method  discussed  previously) , since  these  quantities  are  variables 
within  the  linear  formulation,  and  consequently  are  determined  as 
part  of  the  linear  programming  solution.  An  additional  feature  of 
the  technique  presented  is  that  the  linear  design  constraints  are 
obtained  without  the  use  of  Taylor  series  expansions  or  various 
other  approximations.  In  the  initial  formulation,  some  of  the  design 
constraints  were  nonlinear,  but  by  holding  a small  number  of  parame- 
ters constant  all  constraints  become  linear  in  the  design  variables.  The 
solution  to  the  problem  is  thus  obtained  in  an  iterative  manner,  i.e., 
initial  values  for  the  parameters  are  selected  and  the  linear  problem 
solved;  then  the  parameters  in  the  nonlinear  constraints  are  revised 
and  the  process  repeated.  This  procedure  continues  until  the  parame- 
ters converge  to  their  optimal  values.  Although  the  technique  is 
iterative,  no  series  expansions  like  those  of  [11]  and  [15],  which 
can  be  very  inaccurate,  are  required.  Also,  efficient  linear  program- 
ming methods  have  been  developed  in  this  work  so  that  each  iteration 
of  the  technique  requires  little  computational  time. 


FORMULATION 

The  problem  considered  can  be  stated  as  follows:  given  a framed 
ii  structural  system  with  fixed  configuration  subject  to  externally 

applied  loads,  find  the  design,  i.e.  the  cross-sectional  size  of  the 
♦ members,  such  that  the  weight  of  the  system  is  minimized  subject  to 
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limitations  on  the  stresses  within  the  elements.  In  this  investi- 
gation, the  loading  on  the  structure  was  restricted  to  a single  load- 
ing condition  with  all  forces  and  moments  concentrated  at  the  nodes. 
Several  other  standard  assumptions  were  made  as  well.  All  structural 
systems  were  assumed  to  br-have  in  a linear  elastic  manner  (see  Ref.  12]), 
and  the  stress  in  the  elements  was  taken  as  a linear  combination  of 
both  the  axial  and  bending  stresses  (see  the  Manual  of  Steel  Construc- 
tion [10]  section  1.6).  All  frame  members  were  assumed  to  be  rolled 
wide  flange  steel  shapes,  with  section  properties,  i.e.  cross-sectional 
area,  moment  of  inertia,  and  section  modulus,  approximated  as  linear 
dependent  quantities.  An  approach  similar  to  this  was  used  by  Brown 
and  Ang  [3] . 

The  objective  function  for  the  general  problem  is  thus: 


w = p Z L . A . 
j=l  ^ 


which  must  be  minimized,  where  A.  and  L.  are  the  cross-sectional 

3 3 

t h 

area  and  length  of  the  ' element  respectively,  p is  the  weight 


per  unit  volume,  and  n is  the  total  number  of  elements  in  the  system. 
The  stress  constraints  can  be  written  in  the  form: 

(F.|  |M  I 

5a^  kHa:  i ^ ^=^'2 " 

3 3 

where  a is  the  allowable  stress,  K represents  a constant  of  propor- 
tionality between  the  section  modulus  and  cross-sectional  area,  is 
the  axial  force  in  the  member  and  is  the  maximum  moment  along 

the  length  of  the  member.  Since  it  was  assumed  that  the  external 
forces  and  moments  are  applied  only  at  the  nodes  of  the  structure. 
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then  the  maximum  stress  will  occur  at  one  of  its  ends.  Which  end 
possesses  the  maximum  stress,  however,  is  unknown,  thus  constraints 
are  written  for  both  the  i and  k ends: 

IM, 


iF-jl 
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F.  I 

— J—  + 


OK 


< A. 
D 


j — 1,2, • . • ,n 


(3) 
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o OK 

The  equilibrium  equations  for  the  structure  can  be  obtained  by 
summing  forces  and  moment  at  the  nodes  producing: 


n 

i; 

j=i 


L.  . S . = R. 
Ji  D 1 


i=l,2. 


rin 


(4) 


where  m is  the  total  number  of  free  nodes,  R.  is  a vector  of  exter- 

1 

nally  applied  forces  and  moments  at  node  i,  is  a vector  of  funda- 
mental  member  forces  for  the  j member: 


(5) 
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and  is  a transformation  matrix: 
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where  t.,  = (X.-X.)/L.  and 
11  k 1 1 

ft 

II 

1 
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are  the  coordinates 

of  the 

end  points 

of  the 
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The  compatibility  conditions  can  be  obtained  by  considering  a 


► ♦ 
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general  element  j as  shown  in  Figure  1.  Position  1-2  corresponds 
to  the  initial  position  before  deformation,  while  the  final  position 
l'-2"  is  composed  of  a combination  of  a rigid  body  displacement 
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from  1-2  to  l'-2'  and  a deformation  from  l'-2'  to  l'-2".  The  rela- 
tionship between  local  deformations  and  fundamental  forces  is: 

L . 
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Transforminq  the  local  deformations  to  global  coordinates,  equation 
(7)  can  . written  as: 


C.U.  . 

D.S. 

3 3 

3=1,2 

t • • • 

,n 

(8) 

where  U.  is  a vector  of 

3 

horizontal  and 

vertical 

deformations  and 

rotations  at  the 

ends  of 

the  member , D . 

3 

is 

the  matrix 

defined  in 

equation  (7),  and  is: 
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(9)  I 

0 

0 0 

0 

1 

-1 

The  relation  between  deformations  and  fundamental  member  forces  for 
the  entire  structure  can  be  obtained  by  combining  the  element  equa- 
tions (8)  and  applying  the  boundary  conditions  producing: 

C U = D S (10) 

where  U is  a vector  of  global  displacements  of  the  free  nodes  in 
the  system,  C is  a matrix  of  direction  cosines,  D is  a matrix  con- 
taining the  Dj  matrices  on  the  diagonal  and  zeros  everywhere  else, 
and  S is  a vector  of  all  the  fundamental  forces  in  the  structure. 
The  size  of  C is  3n  by  3m  and  D is  3n  by  3n,  whert.  n is  the  total 
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number  of  elements  in  the  structure  and  m is  the  number  of  free  nodes. 
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The  compatibility  equations  are  formed  by  eliminating  the  dis- 
placements from  the  last  3n  - 3m  equations  of  (10) . This  is  achieved 
by  performing  a Gaussian  elimination  on  the  matrix  C such  that  the 
upper  3m  by  3m  matrix  becomes  upper  triangular  and  the  lower  3n  - 3m 
by  3m  submatrix  contains  all  zeros.  If  these  same  operations  are 
performed  on  the  matrix  D,  then  the  deformations  will  be  eliminated 
from  the  last  3n  - 3m  equations  in  (10) . These  are  the  desired 
compatibility  equations  of  the  form; 

H S = O (11) 

The  matrix  H,  however,  is  a function  of  the  member  properties,  in 
particular,  the  areas  and  moments  of  inertia  which  are  variables 
in  the  design  problem. 

The  initial  formulation  is  now  complete.  Taking  the  fundamental 
member  forces  and  areas  as  variables  the  objective  function  (1)  and 
equilibrium  constraints  (4)  are  linear.  However,  the  stress  constraints 
(3)  are  nonlinear  because  of  the  absolute  values  and  the  compatibility 
equations  (11)  are  nonlinear  since  H is  a function  of  the  areas  and 
moments  of  inertia.  Thus  linear  programming  techniques  can  not  be  used 
on  the  problem  in  its  current  form. 

TRANSFORMATION  TO  LINEAR  FORM 

The  absolute  values  in  the  stress  constraints  (3)  can  easily  be 
eliminated  by  replacing  each  of  the  inequalities  in  (3)  by  four  con- 
straints of  the  form: 
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F . M.  . 

A.  - ^ ^ > 0 

] a Ko  — 

F M . (12) 

A.  + — 1 ^ > 0 

D a KO  - 

F . M . . 

A.  + + ZT^  > O 

j 0 Ko  — 

and  similarly  for  the  end  of  the  element.  The  four  inequalities 

(12)  impose  the  same  stress  restriction  as  does  the  original  one  in- 
volving absolute  values,  but  are  now  linear  in  the  areas  and  funda- 
mental member  forces. 

The  compatibility  constraints  (11)  cannot  easily  be  converted  to 
linear  form;  however,  they  can  be  written  such  that  the  components  of 


the  H matrix  are  linear  functions  of  moment  of  inertia  ratios: 
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j=2 , 3 . . . , n 


(13) 


which  become  the  only  variables  in  this  matrix.  For  example,  the  H 


matrix  for  the  three  member  portal  frame  shown  in  Figure  2 is: 
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where  the  fundamental  member  force  vector  is  ordered  as: 
S = [Fj^  ^2  ^22  ^23  ^3  ^^33  ^34^ 


(14) 


(15) 
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The  constant  K'  is  the  assumed  ratio  of  the  moment  of  inertia 
to  the  cross-sectional  area  for  rolled  wide  flange  steel  shapes. 


The  ratios  6j  are,  of  course,  functions  of  the  areas,  since; 

K'A^  A 

6 . ^ — = — i = — 

3 I.  K'A.  A. 


(16) 


but  if  these  ratios  are  held  fixed,  then  the  compatibility  constraints 
will  be  linear  in  the  fundamental  member  forces.  Thus  the  entire 
formulation  will  be  a linear  programming  problem.  This  suggests  an 
iterative  approach  where  the  ratios  6^  are  initially  selected,  for 
example  all  set  to  one,  and  the  linear  programming  problem  solved 
for  the  areas  and  internal  forces.  Using  these  areas  and  equation 
(16) , the  ratios  are  changed  and  the  linear  problem  solved  again. 

The  process  is  repeated  until  the  solution  of  the  linear  programming 
problem  from  one  iteration  to  the  next  is  approximately  the  same. 

By  solving  several  examples,  it  appears  that  this  technique  does 
converge  to  the  optimal  solution,  and  in  relatively  few  iterations. 


SIMPLEX  METHOD  OF  SOLUTION 


In  order  to  solve  the  linear  programming  problem  at  each  iter- 
ation, the  revised  simplex  method  was  used  [4,8].  This  technique, 
of  course,  requires  all  variables  to  be  greater  than  or  equal  to 
zero,  which  is  not  the  case  for  the  internal  forces  in  the  current 
formulation.  Thus  a transformation  of  variables  was  performed  using 
three  of  the  stress  constraints  for  each  member.  Equivalent  forces 
were  defined  as; 
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which  must  be  greater  than  or  equal  to  zero  in  order  for  the  stress 
constraints  to  be  satisfied.  Equations  (17)  were  solved  for  the  inter 
nal  forces*  F.,  M..,  and  M.,  in  terms  of  the  equivalent  forces  z 
Z... , and  Z.-  and  substituted  into  the  remaining  formulation-  This 

J z ] J 

transformation  has  the  advantage  that  three  of  the  eight  stress 
constraints  for  each  member  are  eliminated,  while  the  number  of 
variables  in  the  problem  remains  the  same. 

Experience  has  shown  that  the  areas  of  the  members  cannoi  t»e 
permitted  to  approach  zero,  since  this  results  in  extremely  large 
moment  of  inertia  ratios  3^  which  produces  instabilities.  Thus  a 
lower  bound  Q was  placed  on  the  areas,  and  a transformation  perfor«i?<i 
so  as  not  to  increase  the  number  of  constraints.  Equivalent  areas 


Aj  were  defined  as: 

A'.  = A.  - Q 
3 3 


j=l , 2 , . . . ,n 


(lit 


Using  equation  (18)  the  areas  were  then  replaced  with  in  the 
formulation,  where  A^  ^ 0.  The  result  is  a linear  programminq 
problem  in  the  equivalent  forces  and  equivalent  areas 

Aj  for  each  member,  which  can  easily  be  solved  using  the  standard 
revised  simplex  method. 

After  solving  the  linear  programming  problem  with  an  initial 


set  of  moment  of  inertia  ratios,  these  ratios  are  changed,  as 
described  previously,  and  a new  linear  problem  is  formed.  Instead 
of  using  the  standard  revised  simplex  method,  with  a phase  I and 
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phase  II,  to  solve  the  new  problem,  a modified  approach  was  used 
which  greatly  improved  efficiency.  The  technique  was  to  assume  that 
the  same  variables  which  were  basic  at  the  optimum  of  the  old  problem 
remained  basic  for  the  new  problem.  First  it  was  necessary  to  calcu- 
late the  inverse  of  the  basis  matrix  B for  the  new  problem.  But  since 
only  a very  few  constraints  changed  from  the  old  problem  to  the  new, 
i.e.  only  the  compatibility  constraints,  and  since  the  inverse  was 
available  from  solving  the  old  problem,  the  new  inverse  was  calculated 
from  the  old  one  using  the  product  form  of  the  inverse.  This  is 
similar  to  the  procedure  used  in  the  revised  simplex  method  itself  in 
changing  the  basis.  The  only  difference  is  that  rows  of  the  matrix 
are  changed  rather  than  columns.  If  a^  is  a vector  corresponding  to 
the  r^^  row  of  the  basis  matrix  which  is  changed  from  the  old  to  the 
new  problem,  then; 

b"^  = b"^  T (19) 

where  B ^ is  the  old  inverse,  B^^  is  the  new  inverse  after  changing 
the  r^^  row  alone,  and  T is  the  matrix: 


where  I is  the  number  of  constraints  in  the  linear  programming  problem. 
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^ V> 

where  the  k variable  is  basic  in  the  r constraint,  and  to 

are  components  of  the  vector  T which  can  be  calculated  from: 

- -T  -1 
T = a B 
r 

Equation  (19)  is  applied  successively  for  each  row  of  the  coefficient 

matrix  corresponding  to  the  compatibility  constraints.  The  result  is 

the  inverse  of  the  basis  matrix  for  the  new  problem. 

,-l  . 


is  determined  the  values  of  the  basic  variables  Xg  can 


Once  B 

be  calculated.  If  all  of  the  variables  are  greater  than  or  equal  to 
zero,  then  the  standard  revised  simplex  method  can  be  used  to  obtain 
the  optimal  solution  of  the  new  problem.  However,  if  one  or  more  of 
these  variables  is  negative  then  a special  pivotal  operation  is  per- 
formed. The  variable  to  leave  the  basis  is  selected  to  be  the  one 
which  is  most  negative,  i.e. 


Xfir  = I ’^Bi 
^Bi"° 


(20) 


The  variable  to  enter  the  basis  is  selected  corresponding  to: 
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— ^ ^ = max 
^rj 
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(21) 


for  c.  - 2.  < O,  a . < 0,  c.  - 2.  > 0 and  a , O,  or  if  this  does  not 

X X i.  X K K rK 

yield  a variable  x^ , then  corresponding  to: 


c . - 2 . 

1 = max 

^rj 


(22) 


for  ^ O and  a^^  < O,  or  if  this  does  not  yield  a variable 


X j , then  finally  corresponding  to: 


c . - 2 . / c.  - 2.  ) 

_J 1 = niin  J i [ 

"rj  1 -ri  / 


(23) 


for  ^ quantities  are  the  elements  of 

the  coefficient  matrix  for  the  linear  programming  problem  and 
are  the  coefficients  of  the  objective  function  as  the  simplex  method 
proceeds.  If  can  be  selected  corresponding  to  the  first  criterion 

(21)  , then  pivoting  will  result  in  x^  becoming  positive,  the  negative 
variable  x^^  being  driven  to  zero,  the  objective  function  decreasing, 
and  all  positive  Cj^  - remaining  positive.  The  second  criterion 

(22)  is  the  normal  approach  used  for  pivoting  in  the  dual  simplex 
method.  It  can  easily  be  shown  that  selecting  a variable  x^  to  enter 
the  basis  is  always  possible  using  one  of  the  three  criteria  (21) , 

(22)  or  (23). 

This  pivotal  procedure  is  repeated  until  all  negative  variables 
have  been  eliminated,  at  which  point  the  normal  revised  simplex 
method  is  used  to  obtain  the  optimal  solution  of  the  new  problem. 
Experience  has  shown  that  this  procedure  is  very  efficient  in  solving 
the  new  problem,  in  most  cases  requiring  only  a few  pivotal  operations. 

EXAMPLES 

The  technique  developed  in  this  study  is  based  upon  the  presump- 
tion that  the  moment  of  inertia  ratios  6j  will  eventually  converge  to 
their  optimum  values.  The  purpose  of  the  first  three  examples  was  to 
show  that  the  linear  approach  could  indeed  determine  the  optimal  design. 

The  results  of  the  linear  method  were  checked  by  solving  the  same 
problems  via  a nonlinear  programming  technique.  The  nonlinear  approach 
consisted  of  the  Fiacco  and  McCormick  penalty  function  (7)  along  with 
the  Hooke  and  Jeeves  direct  search  algorithm  [9] . The  penalty  function 


used  was: 
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2n 


i ' 


p(x,r)  = f(x)  + r I — 


j=l  9. 


(24) 


where  p(x,r)  is  the  penalty  function,  f (x)  is  the  objective  function, 
is  the  inequality  constraint,  and  r is  the  penalty  parameter. 
The  cross-sectional  areas  of  the  members  were  used  as  design 
variables.  The  original  stress  inequalities  (3)  and  lower  bounds  on 
the  areas  (18)  were  the  only  constraints  used.  At  each  step,  the 
fundamental  member  forces  necessary  to  evaluate  the  stress  constraints 
were  calculated  using  an  appropriate  analysis  technique. 

All  design  examples  incorporated  the  follov^ing  constants: 

2 2 

Q = 5.0  in  (32.3  cm  ) lower  bound  on  the  areas 
a = 24  ksi  (165,600  kN/m^)  allowable  stress 
K = 9.0  section  modulus/area 
K'  =75.0  moment  of  inertia/area 

Also  in  the  linear  approach,  the  optimum  design  was  assumed  to  have 
been  obtained  when  the  change  in  the  objective  function  in  two  suc- 
cessive iterations  was  less  than  0.05  percent. 

Shown  in  Figure  2 is  the  first  design  problem  considered,  which 
is  a portal  frame  subjected  to  a set  of  nodal  loads.  The  moment  of 
inertia  ratios  &3  were  initialized  at  1.0  in  the  linear  techni- 

que, and  the  initial  member  areas  for  the  nonlinear  technique  were 

2 

taken  as  90.0  sq.  in.  (580.5  cm  ).  Table  lA  and  IB  show  the  progres- 
sion of  the  linear  and  nonlinear  techniques,  respectively,  towards  the 
optimum  design.  In  the  linear  approach,  we  define  a cycle  as  the  pro- 
cess of  obtaining  the  optimum  of  the  linear  programming  problem  for  a 
given  set  of  moment  of  inertia  ratios,  while  the  minimization  of  the 
penalty  function  for  a particular  value  of  the  penalty  parameter  r is 


considered  a cycle  in  the  Fiacco  and  McCormick  routine.  For  this 
excunple,  the  number  of  cycles  necessary  for  convergence  was  six  for 
the  linear  approach  as  compared  to  seven  for  the  nonlinear  method  of 
solution.  As  can  be  seen  from  Tables  lA  and  IB  the  optimum  values  of  the 
cross-sectional  areas  obtained  from  both  methods  were  virtually  iden- 
tical. The  objective  function,  which  was  taken  as  the  volume  of  the 
structure  rather  than  weight,  was  15,525.0  cu.  in.  (254,610  cm^)  from 

3 

the  linear  approach,  and  15,526.0  cu.  in.  (254,626  cm  ) from  the  non- 
linear solution.  The  total  time  of  execution  was  00:00:36  (hours: 
minutes : seconds)  for  the  linear  technique  versus  00:00:42  for  the 
Fiacco  and  McCormick  routine,  which  indicates  only  a slight  increase 
in  efficiency  for  the  linear  method. 

The  four  member  frame  shown  in  Figure  3 was  used  as  the  second 
example.  Again,  the  starting  point  for  the  linear  technique  consisted 
of  all  the  moment  of  inertia  ratios  being  1.0,  while  the  initial 

2 

areas  for  the  nonlinear  approach  were  taken  as  75.0  sq.  in.  (483.7  cm  ) 
throughout.  Tables  2A  and  2B  reveal  the  route  both  methods  followed  as 
each  proceeded  towards  the  optimum  design.  As  in  the  first  example, 
the  final  design  obtained  from  the  linear  and  nonlinear  techniques  were 
virtually  the  same  with  the  objective  function  being  22,871.0  cu.  in. 
(375,084  cm^)  from  the  linear  solution  as  compared  to  22,892.0  cu.  in. 

3 

(375,429  cm  ) from  the  nonlinear  approach.  For  this  example  the  number 
of  cycles  necessary  for  convergence  was  eighteen  in  the  linear  technique 
versus  eight  in  the  nonlinear  approach.  Although  the  niomber  of  cycles 
was  larger  for  the  linear  solution,  the  total  execution  time  was  much 
less,  with  the  linear  technique  taking  00:00:45  while  the  Fiacco  and 
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McCormick  routine  required  00:03:18.  Thus,  the  linear  programming 
^ approach  exhibited  a substantial  savings  in  computational  effort  ' 

over  the  nonlinear  routine. 

The  last  example  used  for  comparing  the  linear  technique  against 

the  nonlinear  routine  was  the  ten  member  frame  shown  in  Figure  4.  The 

moment  of  inertia  ratios  for  all  the  members  were  again  initialized 

at  1.0  in  the  linear  approach,  while  the  initial  member  areas  for  the 

nonlinear  technique  were  125.0  sq.  in.  (806.3  cm^)  throughout.  Tables  ] 

3A  and  3B  show  the  progression  of  both  methods  towards  the  optimum 

design.  The  final  structures  from  the  two  techniques  are  similar. 

A slightly  lower  value  of  the  objective  function  was  obtained  from 

3 

the  linear  approach,  being  65,647.0  cu.  in.  (1,076,611  cm  ) as 
compared  to  66,191.0  (1,085,532  cm^)  from  the  Fiacco  and  McCormick 
routine.  This  difference  is  due  to  the  fact  that  the  linear  techni- 
que yielded  a design  with  all  members  fully  stressed,  whereas  the 
nonlinear  approach  satisfied  the  convergence  criteria  before  reaching 
the  fully  stressed  condition.  Computational  time  of  the  two  methods 
varied  drastically,  with  the  linear  approach  taking  00:02:57  versus 
01:42:32  for  the  Fiacco  and  McCormick  routine.  This  example  demon- 
strates that  the  nonlinear  approach  is  very  computationaly  ineffic- 
ient when  applied  to  a problem  with  large  numbers  of  design  variables. 

Shown  in  Figure  5 is  a plot  of  computational  time  versus  the 
number  of  members  in  the  system  for  both  of  the  methods  discussed. 

Also  included  in  the  plot  corresponding  to  the  linear  technique  is 
the  computational  time  for  the  15  member  frame  presented  as  the  next 
example.  As  can  be  seen,  the  linear  technique,  compared  to  the 


I 


particular  nonlinear  method  used  in  this  work,  required  much  less 
computational  time  for  larger  structures. 

A fourth  example  was  examined  to  show  the  effect  of  the  initial 
moment  of  inertia  ratios  on  the  final  design.  Shown  in  Figure  6 is 
a 15  member  two  bay  three  story  frame  subjected  to  a single  nodal 
load  condition.  This  frame  was  designed  using  three  unique  starting 
points,  which  will  be  referred  to  as  Cases  I,  II,  and  III.  In  Case  I,  all 
of  the  moment  of  inertia  ratios  were  initalized  as  1.0,  while  Case  II 
started  with  a value  of  5.0  for  all  of  the  member  ratios,  and  Case 
III,  which  used  a bad  starta.ng  point,  was  initiated  with  values  of  the 
ratios  which  differed  greatly  from  the  optimal  values  obtained  in 
the  first  two  cases.  Table  4 shows  the  initial  and  final  moment  of 
inertia  ratios,  final  member  areas,  and  the  values  of  the  objective 
function  obtained  from  the  three  cases.  The  final  designs  are 
almost  identical  with  the  optimum  value  of  the  objective  function 
being  132,805  cu.  in.  (2,178,002  cm^)  , 13.'  , 843  cu.  in.  (2,178,625  cm^)  , 
and  132,846  cu.  in.  (2,178,674  cm^)  for  Cases  I,  II,  and  III  respec- 
tively. The  number  of  design  cycles  and  computational  times  were: 

7 cycles  and  00:08:21,  13  cycles  and  00:11:32,  and  17  cycles  and 
00:13:31  for  the  three  cases.  Thus,  as  expected,  computational  times 
increased  for  starting  points  which  were  further  from  the  optimal 
design,  but  the  increase  is  not  overwhelming. 

The  last  example  considered  was  the  25  member  frame  shown  in 
Figure  7.  Again  two  different  starting  points  were  used,  with  both 
converging  to  approximately  the  same  final  design  as  shown  in  Table  5. 
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It  is  interesting  to  note  that  eight  of  the  member  areas  in  the 

2 

optimal  design  were  at  their  lower  bound  of  5 sq.  in.  (32.3  cm  ). 
Thus,  if  some  or  all  of  these  members  were  eliminated  from  the 
structure  and  the  optimal  design  obtained  again,  it  would  likely 
have  a lower  weight.  The  elimination  of  the  members  might,  how- 
ever, be  influenced  by  other  considerations,  such  as  stability, 
deflections,  etc. 

CONCLUSIONS 

The  linear  technique  presented  in  this  paper  appears  to  be  an 
effective  method  of  optimizing  a framed  structure.  Although  no 
proof  of  convergence  has  been  developed,  the  method  was  able  to 
determine  the  optimum  design  for  all  examples  considered  and  all 
initial  starting  moment  of  inertia  ratios.  It  also  seems  to  be  a 
very  efficient  method.  For  larger  problems,  the  computational  time 
for  the  linear  technique-  .vas  much  less  than  for  the  nonlinear  method. 
Although  other  nonlinear  programming  techniques,  or  improvements 
in  the  one  used,  could  lead  to  decreased  computational  times,  it  is 
felt  from  previous  experience  that  most  nonlinear  methods  could  not 
be  computationally  competitive.  Also,  nonlinear  techniques  have 
other  problems  associated  with  them,  such  as  relative  minim£<,  pre- 
mature termination,  etc.  which  do  not  seem  to  occur  for  the  linear 
method . 

Although  only  simple  loading  conditions  and  stress  constraints 
were  considered  in  this  work,  it  should  be  easy  to  extend  the  method 
to  include  multiple  loading  conditions,  and  non-nodal  loads  (as  long 
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as  finding  the  maximum  stress  in  the  members  does  not  result  in 
complex  nonlinear  constraints) . Likewise  displacement  constraints 
can  easily  be  added.  This  can  be  done  by  continuing  the  Gaussian 
elimination  on  the  force-deflection  equations  (10)  used  to  obtain 
the  compatibility  constraints  (11) . The  goal  would  be  to  reduce  the 
upper  3m  by  3m  matrix  of  C to  the  identity  matrix,  thus  producing 
an  expression  for  the  displacements: 

U = G S (25) 

where  g is  the  upper  3m  by  3n  part  of  the  matrix  D of  equation  (10) 
after  the  Gaussian  elimination  has  been  performed  on  it.  Then  dis- 
placements constraints  of  the  form: 


G S < U 

— max 


(26) 


could  be  included  on  the  problem.  These  can  be  written  such  that 
the  moment  of  inertia  ratios  6^  are  the  only  variables  multiplying 
the  fundamental  forces  S;  therefore,  the  displacement  constraints 
can  be  handled  in  the  same  way  as  the  compatibility  conditions  were 
satisfied  in  method  presented  here. 

With  these  extensions,  it  is  believed  that  the  linear  program- 
ming technique  developed  in  this  research  can  be  an  effective  method 
of  solving  optimal  frame  design  problems. 
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APPENDIX  II.  - NOTATION 

The  following  symbols  are  used  in  this  paper: 
A.  = cross-sectional  area  of  the  member; 


= equivalent  area; 

= elements  of  the  coefficient  matrix; 
= basis  matrix; 


C = matrix  of  direction  cosines; 

Cj  - Zj  = objective  function  coefficients; 

D = force  matrix; 

e.  = unit  vectors; 

1 

Fj  = axial  force  in  jth  member; 

g^  = inequality  constraint; 

H = compatibility  matrix; 

t h 

Ij  = moment  of  inertia  in  j member; 

K = section  modulus/area; 

K’  = moment  of  inertia/area; 

Lj  = length  of  member; 

L. .  = force  transformation  matrix; 

31 

M.  = moment  in  the  jth  member; 

3 

m = number  of  joints  in  the  structure; 

n = number  of  members  in  the  structure; 

p = penalty  function; 
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= lower  bounds  on  the  areas; 

= vector  of  external  loads  at  joint; 

= penalty  function  parameter; 

= vector  of  internal  loads; 

= matrix  used  to  calculate  B inverse; 

= direction  cosines  of  j member; 

= vector  of  nodal  displacements; 

= weight  of  structure; 

= vector  of  variables; 

= equivalent  forces  in  member; 

th 

= moment  of  inertia  ratio  for  j member; 
= vector  used  to  calculate  B inverse; 

= weight  per  unit  volume; 

= allowable  stress; 

= vector  used  to  calculate  B inverse. 


Y 


200*^ 


(0,0)  (100,0) 


Figure  2:  Three  Member  Frame  (1  l<ip  = 4.4  5 kN; 

1 in  = 25.4  mm;  1 kip-in  = 113  N*m) 


Figure  3:  Four  Member  Frame  (1  kip  = 4.45  kN; 

1 in  - 25.4  mm;  1 kip-in  = 113  N-m) 
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: Computational  Time  versus  Number  of  Members 
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Table  2A:  Results  for  Four  Member  Frame  Using 
Linear  Programming  Technique 


Cycles  Moment  of  Inertia  Ratios  Cross  Sectional  Areas  Objective 


Fiacco  and  McCormick 


Member 


Ob j ective 
Function 


Cycle  3 


13.728 
3.324 
4 . 359 
0.994 
0.830 
4 . 028 
1.271 
7 .182 
25.696 


(Sq.  In . ) 


128 .48 
9.36 
38.66 
29.49 
129.20 
154.72 
31.90 
101 . 11 
17.89 
5. 00 


65527.0  in^ 


Cycle  6 

A . 

D 

(Sq.  In.) 

129.71 

15.977 

8 .12 

3.084 

42.06 

3.795 

34  . 18 

0.959 

135.19 

0.867 

149.59 

7.608 

17 .05 

1 . 080 

120.14 

13.723 

9.45 

25.941 

5.00 

65647. 

0 in^ 

Computational  Time  (Hr:Min:Sec)  = 00:02:57] 


2 3 

Note:  1 sq  in  = 6.45  cm  ; 1 cu  in  = 16.4  cm 


Table  3A;  Results  for  Ten  Member  Frame  Using 
Linear  Programming  Technique 


I 


Member 

Starting  Point 

Cycle  3 

Cycle  7 

(Sq.-^In. ) 

A . 

(Sq.^In. ) 

A . 

(Sq.^In. ) 

i 

1 

-- 

125.00 

— 

130.28 

— 

129.48 

2 

1.0 

125.00 

13.226 

9.85 

15.676 

8.26 

3 

1.0 

125.00 

3.083 

42.26 

3.100 

41.77 

4 

1.0 

125.00 

3.382 

38.52 

3.456 

37.46 

5 

1.0 

125.00 

0.965 

135.01 

0.966 

134 . 04 

6 

1.0 

125.00 

0.877 

148.63 

0.878 

147.43 

7 

1.0 

125.00 

6.933 

18.79 

7.758 

16.69 

8 

1 . 0 

125.00 

1.077 

121.00 

1.076 

120.31 

9 

1.0 

125.00 

9.953 

13.09 

9.437 

13.72 

10 

1.0 

125.00 

21.933 

5.94 

25.896 

5.00 

Objective 
Function  , 

— — 

— — — 

67127 

. 0 in^ 

66191. ( 

3 in^ 

Computational  Time  (HriMin: 

Sec)  = 01 

;42:32) 

Note ; 


2 3 

1 sq  in  = 6.45  cm  ; 1 cu  in  = 16.4  cm 


Table  3B:  Results  for  Ten  Member  Frame  Using 
Fiacco  and  McCormick  Routine 


Table  4 : Results  for  15  Member  Frame  Using 
Linear  Programming  Technique 


Member 


Initial 


CASE  I 


Final 


Objective 

Function 


187275.0  in^ 


CASE  II 

Initial  , 

Final 

A . 

3 

Final 
(Sq . In  . ) 

— 

134 .43 

0.907 

148.24 

0.869 

154.64 

4.836 

27.80 

1 . 045 

128 . 63 

26.887 

5.00 

1.037 

129.59 

7.242 

18.56 

0.814 

165.17 

26.887 

5.  00 

1.114 

120.73 

1.210 

111.11 

26.887 

5.00 

1 . 095 

122.72 

26.887 

5.00 

2.512 

53.52 

0.705 

190.72 

26.887 

5.00 

1.117 

120.34 

26.887 

5.00 

1.093 

122.94 

15.314 

8.78 

26.887 

5.  00 

26.887 

5.00 

2.838 

47.37 

0.929 
0.896 
4.869 
1 . 070 
27 .601 
1 . 061 
8 .782 
0.847 
27 .601 


27 .601 
1.131 

27.601 
2.606 
0.720 

27.601 
1 .153 

27.601 


16.035 

27.601 

27.601 

2.829 


2 3 

Note:  1 sq  in  = 6.45  cm  ; 1 cu  in  = 16.4  cm 


138.00 
148.58 
154 . 08 
28.34 
128.93 

5.00 
130.10 
2 


162 . 97 

5.00 

120. 97 
111 . 06 

5.00 

122.00 

5.00 
52.96 
191 .76 

5 . 00 
119.70 

5.00 
123.74 
8.61 

5.00 

5.00 
48.78 


187421.0  in^ 


Table  5:  Results  for  25  Member  Frame  Using 
Linear  Programming  Technique 
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